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Abstract. We have analyzed the properties of dust in the high galactic latitude translucent cloud Lynds 1780 using ISOPHOT 
maps at 100/im and 200/jm and raster scans at 60pm, 80/zm, lOOjum, 120//m, 150/jm and 200/jm. In far-infrared (FIR) emission, 
the cloud has a single core that coincides with the maxima of visual extinction and 200//m optical depth. At the resolution of 
3.0', the maximum visual extinction is 4.0 mag. At the cloud core, the minimum temperature and the maximum lOO/im optical 
depth are 14.9±0.4 K and 2.0 ± 0.2 x 10"', respectively, at the resolution of 1.5'. The cloud mass is estimated to be ISMq. The 
FIR observations, combined with IRAS observations, suggest the presence of different, spatially distinct dust grain populations 
in the cloud: the FIR core region is the realm of the "classical" large grains, whereas the very small grains and the PAHs have 
separate maxima on the Eastern side of the cold core, towards the "tail" of this cometary-shaped cloud. The color ratios indicate 
an overabundance of PAHs and VSGs in L1780. Our FIR observations combined with the optical extinction data indicate an 
increase of the emissivity of the big grain dust component in the cold core, suggesting grain coagulation or some other change 
in the properties of the large grains. Based on our observations, we also address the question, to what extent the 80jum emission 
and even the lOOjum and the 120jum emission contain a contribution from the small-grain component. 

Key words. ISM: clouds - infrared: ISM - ISM: individual objects: L1780 



1. Introduction 

Recent ISM dust models (e.g. Desert, Boulanger & Puget 1990) 
typically have three dust components for which a power law 
size distribution is assumed: the big "classical" silicate or 
carbonaceous dust grains with maximum emission at about 
100-200;um; very small grains (VSGs), which mainly emit at 
A-20- 60/im; and polycyclic aromatic hydrocarbons (PAHs) 
which have emission features in the wavelength range 3-20/im 
and a possible continuum contribution. It has been known since 
the IRAS observations that both the small and the big interstel- 
lar grains contribute to the heating and cooling of the interstel- 
lar medium (ISM), and it is believed that also PAHs contribute 
significantly to the heating process (Bakes & Tielens 1994). 
In a typical translucent cloud, such as L1780, heating takes 
place mainly through photoelectrons that are created when en- 
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ergetic photons from the interstellar radiation field (ISRF) hit 
dust grains. When heated by the solar neighborhood ISRF, 
large grains are able to remain at an equihbrium temperature, 
emitting mostly at far-infrared (FIR) wavelengths. On the other 
hand, when small grains or PAHs absorb an energetic UV pho- 
ton they are temporarily heated to higher temperatures, which 
causes their temperature to fluctuate. While small grains, which 
radiate mostly at shorter wavelengths (below lOOjum), play a 
significant part in both the heating and cooling of clouds, their 
relationship to the large grains, their relative abundance, as well 
as their internal properties, are still incompletely understood. 

CORE and ISO satelUtes extended the infrared observa- 
tional range to longer wavelengths, enabling the measurement 
of the large-grain temperatures and emissivities. Observations 
show that the spectral energy distribution of the interstellar 
dust varies from the diffuse medium to molecular clouds and 
also from cloud to cloud. This probably imphes changes in the 
dust composition, and especially variations in the abundance of 
VSGs and PAHs (Boulanger et al. 1990, Lagache et al. 1998, 
Verier et al. 2000). Variations in the abundances of grain popu- 
lations within a single cloud have also been suggested (Bernard 
et al. 1992, 1993; Rawhngs et al. 2005). Moreover, the ob- 
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served very cold temperatures and changes in the emissivity 
of the large grains also suggest that the properties of the large 
grains change in regions of high dust column density in the in- 
teriors of molecular clouds (Bernard et al. 1999, Stepnik et al. 
2003). In this paper, and in a subsequent paper, where a radia- 
tive transfer model of L1780 will be presented (Ridderstad et 
al., in preparation), we address some of the above questions 
through observations obtained towards the translucent cloud 
Lynds 1780. 

Lynds 1780 (hereafter L1780; also listed as MBM 33 in 
Magnani, Blitz & Mundy (1985) is a high galactic latitude 
(/ = 359.0°, b = 36.7°) cloud with the size of ~ 40' x 30'. 
The brighter central part of the cloud was listed as a separate 
source, L1778, in Lynds' (1962) catalogue. L1780 was classi- 
fied as a cometary globule associated with Loop I by Toth et 
al. (1995). It is located in the dark cloud complex that includes 
the clouds LI34, L169, and LI83 (L134N), and is likely to be 
physically associated with it (Clark & Johnson 1981). We refer 
to these clouds, including L1780, as the L134 complex. The 
distance to the LI 34 complex was determined to be 1 10+10 
pc by Stromgren photometry (Franco 1989) and 140±20 pc by 
Vilnius photometry (Cernis & Straizys 1992). The Nal mea- 
surements by Lallement et al. (2003) indicate a distance of 
~ 100 pc. In this paper, we adopt the distance of 1 10 pc. 

The LI 34 complex is the northernmost extension of the 
Scorpius-Ophiuchus star-forming region and is located at the 
border of the Local Bubble (Kuntz, Snowden & Verter 1997, 
Lallement et al. 2003). In addition to L1780, also other clouds 
in this complex have displaced, cometary-Uke cores. In L183, 
very early phases of star formation are indicated (Ward- 
Thompson et al. 1994, 2000; Lehtinen et al. 2003). No star 
formation has been observed in L1780, although young Ha 
emission line stars with ages a few times 10'' years have been 
observed around L1780 (Martin & Kun 1996). One of these 
Ha stars has the same radial velocity as L1780 (Martin & Kun 
1996), and they may be related to an old SNR expansion shell 
(Toth et al. 1995, Martin & Kun 1996). 

Observations of HI in L1780 have shown that, on its 
Southern half, the cloud has excess emission, which is prob- 
ably caused by the UV-radiation of the bright OB stars in the 
direction of the galactic plane (Mattila & Sandell 1979). Toth 
et al. (1995) discussed in detail the virial equilibrium condi- 
tions and the evolutionary history of LI 780, and suggested that 
the cometary structure of the cloud is produced by the shock 
fronts of supernovae and stellar winds from the Sco-Cen asso- 
ciation of OB stars. They also found that the '^CO core, which 
coincides with the IRAS lOOjum emission core, is in virial equi- 
librium. 

The maximum optical extinction in the cloud has been 
estimated to be 4 mag by Mattila (1986). He also found a good 
correlation between CH column density and optical extinction 
in L1780. The HI and CO distributions in L1780 (Mattila & 
Sandell 1979, Toth et al. 1995) were correlated with different 
IRAS band distributions, which suggested that there are differ- 
ent grain populations present in LI 780, and their distributions 
are distinct from each other and related to the local physical 
conditions. Laureijs (1989) studied L1780 in the optical and 
infrared, and suggested that the optical excess red emission ob- 



served in the cloud by Mattila (1979) is analogous to that ob- 
served in the Red Rectangle (Schmidt, Cohen & Margon al. 
1980). If this interpretation is correct, L1780 is so far the only 
translucent cloud in which this Extended Red Emission (ERE) 
has been observed. 

In this paper, we study the 60jum-200/im emission of L1780 
using ISO observations. We also compare these data with the 
IRAS observations at 12//m and 25jum to investigate the PAH 
abundance within L1780, as the radiation at 12//m is generally 
believed to be due to PAHs, and the 25/im emission is mainly 
from VSGs (Desert, Boulanger & Puget 1990). While the radi- 
ation at 60/vm is due to VSGs, the radiation at ~100/vm-200yum 
can be described as a modified blackbody emission from the 
'classical big' dust grains, which are in equilibrium with the 
surrounding radiation field. It is not entirely clear, to what ex- 
tent the lOOyum emission contains a small-grain contribution; 
we compare our ISO raster scan observations to the current 
ISM dust model to address this question. However, since the 
small-grain contribution at lOO/^m is small, the lOOjum-200/im 
maps and raster scans can be used to effectively trace the do- 
main of the large grains, which are responsible for most of the 
interstellar extinction at optical wavelengths and make up most 
of the dust mass. We use the 100//m and lOOfim ISO observa- 
tions to derive estimates on the dust temperature, column den- 
sity and mass. From the 60/zm, 80jum, 100/vm, llOfim, 150/im 
and 200jum ISO raster scans, the spectral energy distributions 
(SEDs) and temperatures at selected positions of the cloud are 
obtained. 2MASS data are used to derive an optical extinction 
map. The locations of the maxima and the spatial distributions 
of the emission at different wavelengths are compared, and the 
distributions of the three dust components (PAHs, VSGs and 
the large grains) are estimated. We also compare our results 
on LI 780 with the recent observations suggesting that the FIR 
emissivity of dust increases where low dust temperatures are 
observed (Cambresy et al. 2001, del Burgo et al. 2003, Stepnik 
et al. 2003, Kramer et al. 2003). 

2. Observations and data reduction 

2. 1. FIR observations 

The observations were made with the ISOPHOT instrument 
aboard the Infrared Space Observatory (ISO) (Kessler, Steinz 
& Anderegg 1996) satelHte, using the CI 00 and C200 detec- 
tors (Lemke et al. 1996) with the observing template PHT22 in 
raster mode. The data analysis was done with PIA (ISOPHOT 
Interactive Analysis) VIO.O (Gabriel et al. 1997). At the first 
processing level, the detector ramps were corrected for non- 
linearity in the detector response, glitches in ramps were re- 
moved using the two-threshold glitch recognition method, and 
the ramps were fitted with 1st order polynomials. At subse- 
quent levels the signals were deglitched, reset interval correc- 
tion was applied, signals were linearized for the dependence 
of the detector response on illumination, and orbital position- 
dependent dark currents were subtracted. 

Table [J shows the parameters of the individual obser- 
vations. All the maps were calibrated using the FCS (Fine 
Cahbration Source) measurements bracketing the actual mea- 
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Fig. 1. (a) The lOOyum ISO surface brightness map of L1780 at 43.5" resolution. The contour values range from 4.0 to 18.0 
MJy sr"' in steps of 2.0 MJy sr The area used for the background sky value determination and subtraction is marked with 
a rectangle, (b) The 200yum ISO surface brightness map of L1780 at the resolution of 1.5', with contours from 13 MJy sr ' to 
61 MJy sr"' in steps of 6 MJy sr"' . (c) The visual extinction map of L1780 at 3.0' resolution, with contours from 0.8 mag to 4.0 
mag in steps of 0.8 mag, i.e. one step corresponds to about 2cr. The tickmarks point towards decreasing Ay. (d) The 200jum optical 
depth map of L1780 at 1.5' resolution, with contours from 2.0x10""^ to 16.0x10"'' in steps of 2.0x10"''. (e) The temperature map 
of L1780 at the resolution of 1.5', with contours from 15.2 K to 16.8 K in steps of 0.5 K. The tickmarks point towards decreasing 
T. (f) The lOOjum ISO map showing the positions of the IRAS 12/zm (dotted Une), IRAS 25/zm (thick solid line), IRAS 60yum 
(thin solid line), ISO 100//m (dash-dot line) and ISOCAM 6.7jum (large triangle) maxima (Miville-Deschenes 2002), and the 
positions of the detector center of an ISO raster scan (marked with crosses). The cirles show the positions used in calculating 
the averaeed values of the raster scans fsee the text). The black circle indicates the skv nosition used for background subtraction 



M. Ridderstad et al.: Properties of dust in tiie higii-latitude translucent cloud L1780 



5 



Table 1. The parameters of the individual maps and raster scans. Are/is the PHT reference wavelength, Ac is the central wave- 
length, AA is the width, and TDT is the Target Dedicated Time number of the observation. 



Filter A„f A^ AA TDT Raster steps Map size Camera 

fjm fim fim 
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47.1 
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56.9 
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surements. The sizes of the map pixels of the ClOO maps and 
the C200 maps are 43.5" and 89", respectively. A typical sta- 
tistical uncertainty is 0.2 MJy sr ' for a map pixel of the large 
lOOyum map and 1.2 MJy sr"' for a 200//m map pixel. For the 
raster scans, the average uncertainties are 0.4, 0.4, 0.3, 1.5, 0.1 
and 0.2 MJy sr"' for a 60//m, 80//m, lOOjum, 120//m, 150jum 
and 200jL/m map pixel, respectively. The uncertainty in the ab- 
solute calibration is -25 % for the C 1 00 and -20% for the C200 
data (Klaas et al. 2000). 

For the flat-field correction of the large ClOO and C200 
maps a statistical method was applied: the pixel values were 
correlated against the reference pixel at each raster position. 
Rather than comparing the reference pixel only with pixels be- 
longing to the same raster position, a mean of two pixels lo- 
cated symmetrically around the reference pixel was taken. This 
reduces the scatter in the pixel-to-pixel relation caused by sur- 
face brightness gradients. For the flat-fielding of the ISO raster 
scans, a similar procedure was used. The value of the reference 
pixel was calculated as a weighted average using a Gaussian 
that was placed at the position of the studied pixel and had a 
FWHM of 2.0'. For both the large ISO maps and the raster 
scans the reference pixels were chosen to be the pixels number 
8 and 3 in the ClOO and C200 rasters, respectively. The data 
from pixel no 3 of the ISO 60jum raster scan was of such bad 
quality that it was replaced by the average of the values of the 
surrounding raster pixels. 

It is known that the primary intensity calibrators, the fine 
calibration source (FCS) measurements used to derive the re- 
sponsivities of the detectors, suffer from signal transients in the 
case of the ClOO camera. Lehtinen et al. (2001) showed that 
even after these transients were corrected using the PIA sig- 
nal drift interface, there remained a difference of 20% between 
the lOOyum ISO and DIRBE surface brigtnesses. Therefore, the 
ClOO measurements on L1780 were compared with DIRBE. 

IRAS data was first scaled to DIRBE scale using the rela- 
tion /DiRBE(lOOjum) - 0.73 X /iRAs(100//m) that was derived 
for a 6° diameter circular region around LI 780. A linear fit 
ISO vs. scaled IRAS values gave /iRAS(DiRBE)(100jum) = 0.61 x 
^iso(lOOjum). In the comparison both DIRBE and ISO obser- 
vations were color corrected assuming a spectrum ByiTdu^i = 
17 K)v^. The larger lOOjum map was scaled to the DIRBE sur- 



face brightness scale using this relation, and the 100/vm raster 
scan was adjusted to the same level with the scaled raster map. 

The 6Qfj.m raster scan values were similarly scaled to 
DIRBE via IRAS, using the relations /DiRBE(60jum) = 0.69 x 
/iRAs(60yL/m) and /iRAS(DiRBE)(60yum) = 0.86 x /iso(60yL/m). No 
color coiTection was applied since at 60yum the dust spectrum 
is rather close to a flat spectrum, v/y = constant. We also com- 
pared the values calculated from the zodiacal light model of 
Good (1994), to the corresponding empty-sky values of the 
60//m raster map (raster position within the circle marked in 
Fig- Ed that were color corrected for a 270 K black body spec- 
trum. This gave the relation /„,ojf/(60jum) = 0.89 x /iso(60//m). 
We took the average of the two scaling factors, obtaining 
^rrae(60/vm) = 0.88 X /iso(60//m) as the final scaling relation. 

The 80//m data was scaled by linearly interpolating the 
scaling coefficients for the 60jum and lOO/^m ISO data, which 
gave 0.71 as the scaling coefficient for 80//m. 

The IRAS 12/im and 25yum Infrared Sky Survey Atlas 
(ISSA) map data were re-scaled to DIRBE scale using the rela- 
tions /DiRBE(12//m) = 1.06 X /iRAs(12//m) and /DiRBE(25jum) = 
1.01 X /iRAs(25yum), as given in the ISSA Explanatory 
Supplement (Wheelock et al. 1994). 

The ISO 200yum map surface brightness values were also 
compared with DIRBE data. There are only 4 DIRBE pixels 
within the 200jum map area. The DIRBE surface brightness 
values at 100//m, 140yum and 240/vm were fitted with a mod- 
ified blackbody function with A^^ emissivity law in order to get 
interpolated values at 200/vm. The ISO 200jum map was con- 
volved with a scan-averaged DIRBE beam and the ISO surface 
brightness values were color corrected using the temperatures 
derived from fitting the DIRBE values. A linear fit forced to go 
through the origin gives /DiRBE(200//m) - 0.98 x /iso(200yum), 
which is well within the estimated uncertainties of the ISO data. 
Since this fit is made based on four points only, which results in 
considerable uncertainty, and the C200 detector does not suffer 
from such problems as described above for the ClOO detector, 
we retain the FCS based calibration for the 200//m and other 
C200 ISO data (the 120yum and 150yum data). 

The sky area used for subtracting the background values in 
the 100/vm and 200yum maps is shown in Fig.^. The back- 
ground sky values for background subtraction for all the six 
raster scans and the IRAS data used were determined at the 
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Fig. 2. The lOOyum vs. the 200;um ISO map values. The solid 
lines show the surface brightness ratios at the indicated tem- 
peratures, assuming an/3 - 2 emissivity power law index. The 
points corresponding to the sky area of the ISO lOOyum surface 
brightness maximum (6.25' x4.0', corresponding to the highest 
level contour in Fig.[^) are marked with squares. 



same position outside the cloud (Fig.^); this way we obtained 
a common zero level for comparison with the raster scans. To 
compare the six raster scans with each other, we calculated 
average values at 16 evenly-spaced positions along the CI 00 
scans, and at 13 positions along the (shorter) C200 scans (see 
Fig. (2" for these positions; three ClOO positions in the west fall 
outside the figure). The averaging was done using a Gaussian 
weight function with FWHM=4.5'. This resolution also corre- 
sponds to the resolution of the ISSA maps. 

2.2. Extinction data 

An extinction map of LI 780 was derived using the J, H and 
band magnitudes obtained from the 2MASS archive. The 
extinction measurements are based on the near-infrared (NIR) 
color excesses of stars visible through the cloud. We applied the 
optimized multi-band NICER technique of Lombardi & Alves 
(2001), which is a generalization of the traditional color excess 
method (using data on two bands only), to derive the NIR color 
excesses. 

Stellar density of the detected stars at all three bands is 
fairly constant over the map, with an average of about one 
star per square arc minute. As the reference field, represent- 
ing an area without significant extinction, we selected a 1° di- 
ameter circular area located close to L1780 and at the same 
galactic latitude. The coordinates of the centre of this area are 
R.A.(J2000)= 15''22™00^ Dec.(J2000)= -6°05'00". The col- 
ors of the stars in the reference field have the mean values and 
standard deviations of (7 - //)„ = 0.49 + 0.19 and {H - K)o = 
0.12 + 0.23. As the effective wavelengths of the 2MASS J, 
H and bands we use 1.25;um, 1.65//m and 2.17/im, respec- 
tively (Kleinmann et al. 1994). For the ratios of visual extiction 
to color excess we used the values Av/E{J-H) = 8.86 and Ay/ 



E(H - K^) - 15.98, which correspond to an extinction curve 
with Rv = 3.1 (Mathis 1990). 

The extinction value in each map pixel was derived from 
the individual extinction values of the stars by applying the 
sigma-clipping technique of Lombardi & Alves (2001). In or- 
der to obtain a sufficient signal-to-noise ratio in the Ay map, the 
individual extinction values were averaged using a Gaussian 
with FWHM=3.0'. 

The extinction map has two eiTor sources, the variance of 
the intrinsic colors {J-H)q and {H-K)q, and the variance of the 
observed magnitudes of the field stars. The former dominates 
with a Icr error of ~ 0.42mag per pixel, while the latter gives 
a Icr error of ~ 0.23mag per pixel, resulting in a typical total 
error of 0.47 mag in Ay. 

The map of visual extinction is shown in Fig. QJ. At the 
resolution of 3.0' the value of maximum extinction is 4.0 mag. 

When the extinction data is compared with the lOO/^m and 
IQQjjm surface brightness, temperature and optical depth maps, 
the latter are convolved to this resolution. 

3. Results 

3. 1. General structure of tlie FIR emission from 
L1780 

The 100//m and IQQfim surface brightness maps are shown in 
Fig. and Fig. In the far-infrared, L1780 has a clearly 
visible brighter central core in the West, which coincides with 
the region of maximum optical extinction. On the opposite side, 
the cloud has an extended, tail-like structure, and an overall 
cometary shape. 

The relation between the 100//m and the 200//m emission 
is shown in Fig.|2] The lOOyum emission is well correlated with 
the lOOjjm emission, and the relation is linear up to the high- 
est surface brightness values, where there is a slight flattening 
trend. The solid lines in the figure show the surface bright- 
ness ratios at the indicated temperatures, assuming a modi- 
fied blackbody function with dust emissivity power law index 
yS = 2. 

As can be seen in Fig. the IRAS 60/im maximum is 
distinct from the 100;um/200;um maximum, and is situated East 
of it. The 25/im, the 12//m and the 6.7jum maxima are located 
still further to the East, reaching towards the SE edge of the 
cloud. The 100//m and IQQfim maps do not show any structure 
at the locations of these maxima. 

The '^CO core of L1780 is in virial equiUbrium (Toth et al. 
1995). However, no point sources, neither those from the IRAS 
Point Source Catalog nor any new ones, are detectable in the 
ISO maps. Also, a J-H vs. H-K color excess plot made using 
2MASS data reveals no stars with color excess >1.0 mag in 
LI 780. This leads to the conclusion that there is no (detectable) 
star formation activity in LI 780. 

3.2. Surface brightness ratios from ISO raster 
scans and IRAS observations 

In Fig. [J, the positions of the ISO raster scans, which go 
through the cloud from West to East and back, are shown. 
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Fig. |3p shows the intensities of the raster scans averaged (us- 
ing Gaussian beams with FWHM=4.5') at the positions marked 
with circles in Fig-^^. 

All the raster scans go through the high density core of 
LI 780 and the regions of the maximum llfim, 25//m and 60jum 
emission. The maxima of the 120/im-150//m emission coincide 
with the maximum of the 200//m emission, whereas the posi- 
tions of the 60/im and the 80//m maxima are shifted towards 
East. The shape of the lOOjum brightness profile is between the 
shapes of the 80/im and the 120jum, especially in the region of 
maximum emission. 

Fig.|3t also shows the values of the IRAS 12jum and 25/vm 
observations calculated at the same positions as the raster 
scans. The maximum of the IRAS 25yum emission is located 
between the 60jum maximum and the 12/im maximum on the 
Eastern side of L1780. 

Fig-Eb shows the intensity ratios 12jum/60jt/m; /y/lOOyum, 
where ly is 12yum, 25jum or 60yum; and /v/200jum, where ly 
is 60//, 80//m, 100//m, 120//m or 150//m, plotted through the 
cloud from East to West (the intensities have been calculated 
as in Fig.|3^). The cold core of L1780 shows up as a minimum 
in all the curves. The ratios /y/lOOjum and 60//m/200yum rise 
steeply towards the Eastern side of the cloud, indicating an in- 
creased presence of the small-grain dust component. The ratio 
12/im/60//m peaks at the maximum position of the 12jum emis- 
sion in LI 780, indicating the presence of R\Hs in this loca- 
tion. The ratios 80jum/200//m and 100//m/200jt/m show a deep 
minimum at the core of 1780, while the ratios 120jum/200//m 
and 150/im/200/Ym show only a slight drop at this location, the 
150yum/200//m ratio being almost constant. 

In Fig. |4] the 60-150jum surface brightness values of the 
raster scans are plotted against the 200//m scan values. All the 
wavelengths coiTelate well with the 200/vm intensity. However, 
with the exception of the ISOjum vs. lOQ/um relation, all curves 
clearly show signs of flattening close to the 200jum maximum. 

3.3. Dust temperature and the spectral energy 
distribution 

The temperature of dust has been derived, pixel by pixel, from 
the lOOyum and 200jum ISO surface brightness maps by using a 
modified blackbody function 

1(A) oc A-f^BiA, Tdust), (1) 

where /3 is the dust emissivity power law index, and B(A, Tdust) 
is the Planck function. We assume fi - 2. For the temperature 
calculation, the lOOjum map was convolved to the resolution 
of the 200//m map. The surface brightness values were itera- 
tively color corrected, using the color correction coefficients 
based on previous temperature determination, until the differ- 
ences between two consecutive iterations were below 0.1 K. 
The temperature map is shown in Fig.[^. At the resolution of 
1.5', the minimum temperature is 14.9+0.4 K. 

The temperature of LI 780 has also been determined using 
the ISO raster scans. Fig. |5l shows the SEDs at three positions 
along the scan: Figs.|5t,|5]3 andlsj; correspond to the locations 
of the IRAS ll^im, the IRAS 60yum and the ISO lOOyum emis- 
sion maxima, respectively (see Fig. for the locations of these 



maxima in L1780). In addition to the ISO raster scan values at 
60-200/vm, the IRAS 12//m and 25yum values are plotted in the 
figure. Modified blackbody curves with dust emissivity index 
P - 2 have been fitted to the four longest wavelengths, result- 
ing in temperatures 15.0+0.4 K, 15.0±0.4 K and 14.1+0.4 K, 
for Figs. [5^, Isj) and|5j;, respectively. These temperatures are 
in agreement with the values obtained using the large 100/im 
and 200yum ISO maps. In the fitting, the statistical errors (from 
PIA) plus filter-to-filter errors of 10% have been used for the 
ISO C200 values (120//m, 150//m and 200//m). The error used 
for the 100/vm values is the statistical error plus a 10% error 
resulting from the uncertainty of the scaling of the 100//m to 
DIRBE. The average value of all the fits was slightly below 
2. In Fig.|5]the 80//m are ~2cr above the fitted curves. The ratio 
of the mid-infrared intensity (average of the 12//m and 25yum 
intensities) to the FIR intensity (maximum of the blackbody fit 
at the three longest wavelengths) is 0.045 for the IRAS \2^m 
maximum region (Fig.|5^), 0.029 for the IRAS 60jum maximum 
(Fig.Eb) and 0.013 for the ISO 100//m maximum (Fig.El;). All 
the surface brightnesses and errors used in Fig. [S] are given in 
TableEl 

3.4. FIR-optical deptti and cloud mass 

The optical depth at 20Qfj.m has been calculated from the equa- 
tion 

/(200//m) 



r(200A/m) = 



(2) 



B(200A/m,rdust)' 
which is valid for optically thin emission and where / is the 
observed intensity and B(rdust) the intensity of a blackbody 
at the temperature Tdust- The optical depth map is shown in 
Fig. QJl. The maximum value of the 200/vm optical depth is 
2.0+0.3 X 10"^ at the resolution of 1.5'. 

To calculate the mass of the cloud we first derived the ratio 
between the FIR optical depth and the hydrogen column den- 
sity, i.e. the average absorption cross section per H-nucleus 

T(A) 



cr"(A): 



NiH) 



(3) 



The NIR extinction can be used to estimate the total hydrogen 
column density, N(H+H2) in L1780. We adopt the value N(H+ 
H2)/E(B - y) = 5.8 X lO^icm^-mag-i, valid for diffuse clouds 
(BoMin, Savage & Drake 1978), together with Av/E(B -V) = 
3.1 to obtain N(H + //2)/Ay = 1.87 x lO^'cm-^mag"'. The 
relation between Ay and T(200;um) is shown in Fig.|6t. A linear 
fit, T(200fim) = 3.0 X 10""* X Ay - 6.4 x 10"^ also shown in the 
figure, gives cr^(200;um) = 1.4 x 10"^^ cm^ per H-atom for the 
absorption cross section, which is in agreement with values for 
diffuse ISM by Dwek et al. (1997) (1 .4x lO^^^cm^ per H-atom). 
The relation in Fig.|6^ deviates from linearity in the region of 
highest density in L1780; see Section 4.3 for discussion. 

The total mass (gas plus dust) of L1780 has been calculated 
using the equation 

T(200jum) 



M, 



FIR 



-D othjU, 



(4) 



crff(200//m) 

where r2oo/im is the 200//m optical depth, cr''(200jum) is the 
absorption cross section per H-nucleus for which we have used 
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Table 2. The surface brightness values and corresponding eiTors used in Fig.|5l The three raster positions are given according to 
Fig.|3] the raster positions no. 4, 6 and 9 correspond to the ll/im, 60jum and 100//m emission maxima of L1780, respectively. 





Raster position 4 


Raster position 6 


Raster position 9 


A 


/(v)[MJy/sr] 


(5/(y)[MJy/sr] 


/(y) 


6I(v) 


i(y) 


61(v) 


12^ni 


1.06 


0.21 


0.85 


0.17 


0.52 


0.10 


25/im 


1.55 


0.31 


1.36 


0.27 


0.95 


0.19 


60/jm 


2.56 


0.65 


2.92 


0.69 


2.24 


0.62 


80pm 


4.48 


0.87 


6.09 


1.06 


5.94 


1.05 


lOOyum 


7.27 


1.04 


9.78 


1.34 


10.69 


1.44 


120pm 


17.76 


1.93 


23.39 


2.51 


28.86 


3.09 


150pm 


26.90 


2.84 


35.19 


3.69 


46.84 


4.90 


200pm 


24.88 


2.73 
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Fig. 7. (a) ISO 200//m and (b) lOOyum surface brightness values 
versus Ay. 



the value derived above, D is the distance, niu is the hydrogen 
mass, and p is the mean molecular weight. The mass has been 
derived by summing up the pixels of the T(200jum) map. The 
mass thus obtained for L1780 is ISM©. 



In the Ay map, there seems to be another region of higher 
visual extinction, an arc-like structure, on the Eastern side of 
the 100pm core. If real, this feature, together with the cometary 
shape of LI 780, could result from the propagation of a shock- 
wave from the SW direction through the cloud as has been sug- 
gested by Toth et al. (1995). In fact, the arc-like feature could 
also be a part of a ring-shaped structure, in which case an in- 
teresting analog of a possible dense-core formation would be 
found in the Globule 2 of the Coalsack cloud (Lada et al. 2004). 
However, the extinction values may be biased because of vari- 
ations in the stellar density. This hampers the detection of mor- 
phological features in the Ay map at this resolution, which 
is close to the detection limit {AAy ~ 2cr in the contours in 
Fig.Ct). 



Figs.|6]and|8]show the relations between Ay, T(200pm) and 
T^dust- The values from the locations of the IRAS 12pm, IRAS 
60pm and ISO 100pm maxima are marked with triangles, dia- 
monds and squares, respectively. The sizes of sky areas of the 
maxima have been chosen according to the highest-level con- 
tours visible in Fig.^, and are 12.5' x 3.0', 11.3' x 4.0', and 
6.3' X 4.0' for the 12pm, 60pm and 100pm maximum, respec- 
tively. The error bars in the figures show typical errors. The er- 
rors of the temperature and the 200pm optical depth get larger 
for the points near to the borders of the maps, since near the 
borders of the ISO FIR maps the background-subtracted sur- 
face brigthness values, used for calculating Tdust and T(200pm), 
are more uncertain. For the plots in Figs. |6l and |8] most of the 
near-border points have been removed. 



3.5. Comparison of FIR and extinction data 

At large scale, the far-infrared emission in L1780 coiTelates 
well with the optical extinction. The 100pm, 200pm and Ay 
maps in Figs. [^, [Q) and Q]Z! show that the 100pm and the 
200pm emission maxima coincide with the location of the 
highest visual extinction. The relations between the 100pm and 
the 200pm emission and the visual extinction are linear (Fig.0. 



Fig. nil shows the temperature decrease in the densest part 
of the cloud: as the optical depth increases, the temperature 
decreases from ~15.8 K in the regions of the 12pm emission 
maximum (PAHs) and the 60pm maximum (VSGs) to ~14.9 K 
in the 100pm maximum region (big grains). In Fig.|6]3, simi- 
lar dependence is shown for the temperature as the function of 
the visual extinction. Fig.|6^ shows increased emissivity in the 
regions of high Ay in L1780 (see Section 4.3 for discussion). 
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Fig. 6. (a) The 200/^m optical depth versus the visual extinction Ay, at 3.0' resolution. For the linear fit, see Section 3.4. (b) 
Temperature versus the visual extinction, (c) T(200jum)/Av versus Ay. (d) Temperature versus the 200jum optical depth, at 1.5' 
resolution. For (a), (b), (c) and (d), squares, diamonds and triangles mark the points corresponding to the ISO lOOjum, IRAS 
60/im and IRAS 12//m surface brightness maxima in the sky, respectively. The error bars show typical errors. 



Table 3. The average and maximum values of the infrared colors along the ISO raster scans in L1780. The solar neighborhood 
(SN) values are from Boulanger et al. (1990). 



Color 


L1780avg 


L1780max 


L1780min 


SN 


12/100 


0.10 


0.20 


0.05 


0.042 


25/100 


0.16 


0.26 


0.07 


0.054 


60/100 


0.31 


0.42 


0.21 


0.21 


12/25 


0.59 


0.85 


0.29 
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Fig. 8. T(200jum)/Av versus T. The error bars plotted in the up- 
per right corner show the magnitudes of typical errors. The 
curves show the emissivity increment to eight-fold (dashed 
line), four-fold (dash-dot line) and no increment (dotted line), 
when the relative amount of the cold component increases from 
0.2 to 1 (from right to left), according to the model by del Burgo 
et al. (2003). 



4. Discussion 

4. 1. Different dust populations from FIR 
observations 

The Southern side of LI 780 is facing both the Galactic plane 
in the South and the Upper Scorpius (USco) association of OB 
stars in the SW direction, and is, therefore, subject to more in- 
tense ISRF. However, the variations in the emission of L1780 
at different wavelengths are more prominent in the E-W direc- 
tion than on the S-N axis. Fig. [0^ shows that the 12/im emission 
maximum (indicating the presence of PAHs) is located on the 
SE side, in the "tail" of L1780. The 60//m (VSGs) emission 
maximum is located to the West from the 12jum maximum, 
and finally, the lOOjum maximum (big grain emission) is on 
the Western side of the cloud. 

In Fig. [3^, the 60yum emission curve clearly rises towards 
East and has its maximum on the Eastern side of the 200/vm 
maximum. Fig. |3j5 shows that the 60//m/100jum ratio rises 
rather steeply from its minimum towards the East side of 
L1780. In Table|3 it can be seen that the minimum value of this 
color ratio (reached in the cold core of L1780) is the same as 
the solar neighborhood (SN) value by Boulanger et al. (1990). 
Laureijs, Clark & Prusti (1991), who made IRAS observations 
on the LI 34 complex, indicate in thek study of LI 780 that the 
60/im emission diminishes in a narrow transition layer around 
the 100//m maximum. Our data do not show any abrupt change, 
but the decrease in the 60jum/100//m ratio is clear: the ratio is 



Fig. 9. Comparison between the 60jum ISO raster scan and the 
80-120jum raster scans from which the 200/im emission (the 
large-grain contribution) has been subtracted. 

halved in the cold core region relative to the Eastern edge of 
the cloud. It is not likely that this reduction would be due to 
attenuation effects alone: in the model by Bernard et al. (1992) 
the effect of radiation attenuation on the ratio 60um/100yum as 
Av increases from 1 to 4 mag (between different ID models) is 
only half of the decrease of the ratio observed in LI 780. 

Laureijs, Clark & Prusti (1991) attribute the change in 
60um/100jum to the formation of icy mantles on grains. Stepnik 
et al. (2003) shows that the depletion of VSGs through grain- 
grain coagulation can reduce the 60um/100/vm ratio, and, in 
addition, cause an increased FIR emissivity of grains, which 
is also observed in L1780. Fig. |6l; shows the increased val- 
ues of T2oo/Av in the region of the lOOjum emission maximum 
(marked with squares in Figs.|S]3 and|5J;). On the other hand, 
if taken to be traced by the T20[)/Av ratio, no difference in the 
overall size distribution of the large grains between the 60/jm 
and 12jum maxima is visible in Fig. |6l;, although the 6Qfj.m 
maximum is located deeper in the cloud. 

Figs. |3^ and (Sj) show that the emission from VSGs and 
PAHs is concentrated on the Eastern side of LI 780, and that 
the PAH emission clearly is distinct from the VSG emission. 
Also the MIR/FIR ratios from Fig. |5] indicate increase of the 
12jL(m-60/zm emission toward the East. For the 12jum/100jum 
ratio, comparison with the results of the model of Bernard et al. 
(1992) suggests that radiative transfer effects account for less 
than one third of the decrease of the ratio in the LI 780 core. 
Table|3]shows that the average colors 12//m/100yum and 25jum/ 
100/vm are more than twice the SN values. At the scan posi- 
tions of the 12fj.m and 25yum maximum emission, the values of 
12jum/100/vm and 25//m/100yum are 0.15 and 0.18, correspond- 
ing to 3.6 and 3.3 times the SN value. The maximum ll/jm/ 
lOOjum and 25/im/lOO/im ratios observed in L1780 are both 4.8 
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times the SN value. These values can be compared with those 
of two clouds from Chameleon and p-Ophiuchi, which are 5.5 
and 3.3 times the SN value for PAHs (using the I2fim/ IQOjjm 
ratio), respectively, and 2.3 and 2.3 times the SN value for 
VSGs (using the 25yum/100//m ratio), respectively, (Bernard, 
Boulanger & Puget 1993). The color ratios observed in L1780 
indicate that there is an overabundance of PAHs and VSGs on 
the Eastern side of the cloud. Further support to our observa- 
tions on PAH abundance is provided by the strong decrease 
of the 6.5/100//m ratio in the FIR core of L1780 observed by 
Miville-Deschenes (2002) using ISOCAM and IRAS observa- 
tions. 

We conclude that, although attenuation effects may be in- 
volved, they are not sufficient to explain the observed color ra- 
tios, which indicate true variations in the abundances of dust 
grains from different populations. To estimate the effect of ra- 
diation attenuation on the observed color variations in L1780, 
model calculations are needed (Ridderstad et al., in prepara- 
tion). 

4.2. VSG contribution at ttie wavelengttis above 
SO/um 

Radiative transfer modelling performed using the current ISM 
models (including PAHs, VSGs and big grains) indicates that 
the emission from VSGs is still significant above 60yum, and 
contributes from 15% to 23% of the total emission at 100 
fim when Ay ranges from 1 to 4 mag (Desert, Boulanger & 
Puget 1990). Laureijs et al. (1996) made ISOPHOT observa- 
tions at 60jum, 90jum, 135//m and 200yum on a small cloud in 
Chameleon, which, like L1780, has a reduced 60yum/100/im ra- 
tio in the centre, and indicated that the 90/vm emission has an 
excess contribution from dust grains emitting at 60fj.m. 

In Fig.|9] the 200jum raster scan has been assumed to con- 
tain emission solely from the big-grain dust component, which 
then has been subtracted from the 80jum, lOOjum and 120/vm 
raster scans after scaling the latter using a ratio /^^//foo/im deter- 
mined at the cloud center, i.e. at the 200/zm maximum position 
(position 9 in Figs.|3t): 

My = /v - — X 7200^™ MJy/sr. (5) 

^200//m 

The remaining part, A/y, represents excess that can be at- 
tributed to small grain emission, assuming that no temperature 
change of big grains resulting from the increase of A v is taken 
into account. In reality, the A/y always contain also the effect 
from radiation attenuation; thus, Fig.|9]gives an upper hmit to 
the VSG contribution. Fig.|5]shows that the 200yum-subtracted 
intensities peak at the maximum position of the 60/im emis- 
sion. Our results indicate that, at the maximum position of the 
60//m emission, at most 12% of the 120jL/m emission, 22% of 
the 100/^m emission and 28% of the 80/im emission is from 
VSGs. When it is taken into account that the ISO filter band- 
passes overlap, this is consistent with the current ISM dust 
model (see Desert, Boulanger & Puget 1990). 



4.3. Increased FIR emissivity in tlie dense core of 
L1780 

Many recent studies indicate a change in the emissivity of inter- 
stellar grains in dense molecular regions (Cambresy et al. 2001; 
del Burgo et al. 2003; Dupac et al. 2003; Kramer et al. 2003; 
Stepnik et al. 2003; Cambresy, Jarrett & Beichman 2005). Also 
in our study, Fig.|6^ shows that in the cold core of L1780 (the 
region of the FIR emission maximum and minimum temper- 
ature Tdust) above Ay x 3.5, the relation between T(200/im) 
and Ay deviates from linearity, indicating an increased emis- 
sivity of the big grains. This increase in the ratio T(200jum)Mv 
in the 100/im emission maximum area is also visible in Fig.|6j;- 
The emissivity observed in the cold core of L1780 in Fig. [6^ 
is ~1.5 times the values observed elsewhere in the cloud. This 
value is in agreement with the results of Cambresy, Jarrett & 
Beichman (2005), who considered the whole galactic anticen- 
ter hemisphere and found that, for regions with Ay > I mag, 
the ratio of FIR optical depth to NIR extinction is Ay (FIR)/ 
Ay(gal)=l. 3 1+0.06. 

Del Burgo et al. (2003) found FIR emissivity changes in 
their 60-200/im observations of eight regions mostly belonging 
to quiescent high latitude clouds with Ay ~ I - 6 mag. Since 
LI 780 is also a translucent cloud, their results are of particular 
interest here. Del Burgo et al. (2003) presented a model with 
two big-grain dust components: a warm component at T=17.5 
K and a cold component at T=13.5 K. A coefficient e gives 
the factor by which the emissivity of the cold component is in- 
creased. In Fig. |8] the curves showing the emissivity change 
according to the model by del Burgo et al. (2003) are shown 
for comparison with our data. Comparison with the results of 
del Burgo et al. (2003) indicates that the emissivity in the core 
of L1780 is about twice of that for the diffuse ISM, which is in 
agreement with the value derived from the relation in Fig. [6^. 
Although the uncertainty in Fig. |8] increases with increasing 
temperature, it is clear from the average errors (shown by error 
bars in the figure) that there is an increasing trend towards high 
Tioo/Ay values. It must be noted that in this kind of plot, the 
dependance of the optical depth on the dust temperature (since 
the former is calculated using the latter) can cause false corre- 
lation similar to the supposed emissivity increase. However, the 
errors in the plot in Fig.|8]diminish towards high optical depth 
and low temperature. We also have the advantage that the mea- 
surements are from a single cloud, which eliminates the effects 
resulting from calibration eiTors between separately calibrated 
maps. 

The ratio T2ooMv increases with grain size and porosity 
(Cambresy et al. 2001), and the change in FIR/submm emissiv- 
ity has been attributed to the coagulation of grains into larger, 
fluffy particles in dense, cold molecular regions (Cambresy et 
al. 2001, Stepnik et al. 2003). It is interesting that an inde- 
pendent indicator, the decrease of the 60yum/100;um ratio to- 
wards the dense core of many clouds, including L1780, has 
also been attributed to grain coagulation (Bernard et al. 1999, 
Stepnik et al. 2003). While the emissivity change in L1780 is 
seen above ~3.5 mag, Stepnik et al. (2003) derived the thresh- 
old value Ay = 2.1 + 0.5 for the possibly coagulation-induced 
changes (reduced 60//m/100jum ratio and increased emissivity) 
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in the properties of the grains to occur. They noted that these 
changes in the properties of dust grains may be part of a general 
transformation process that dust undergoes in the dense ISM. 
Moreover, Cambresy, Jarrett & Beichman (2005) indicate that 
large, fluffy grains may be conomon all over the galaxy wher- 
ever Ay > 1. 

5. Conclusions 

The analysis of the ISO lOO/im and 200jum maps and 60jum, 
80//m, lOOjum, 120jum, 150//m and 200//m raster scans of 
LI 780, combined with the IRAS 12//m and 25yum data and the 
visual extinction map (based on NIR J, H and band color 
excess data), leads to the following conclusions: 

- The lOOyum emission is well correlated with the 200/zm 
emission throughout the cloud. The cloud has a single core, 
revealed by the lOOyum and 200;um emission, the FIR opti- 
cal depth, and the visual extinction. 

- The spatial distributions of the 12/im, 25jum and 60/im 
emission differ significantly from the emission at longer 
wavelengths. This indicates the presence of separate dust 
components with different physical properties and spatial 
distributions. 

- The maximum values of the color ratios 12yum/100/zm, 
25jum/100jum and 60//m/100/im, which are 4.8, 4.8 and 2.0 
times the solar neighborhood values, indicate an overabun- 
dance of PAHs and VSGs in LI 780. 

- The Ay map gives a maximum visual extinction of 4 mag 
at 3.0' resolution. The visual extinction is well correlated 
with the FIR emission at large scale. 

- The cold core of L1780 has the minimum temperature of 
14.9+0.4 K and the maximum 200yum optical depth of 
2.0+0.3 X 10"^ at 1.5' resolution. 

- The value of the absorption cross section per H-atom at 
200;um has been estimated to be cr'^(200;um) = 1.4 x 
lO^^^cm^ per hydrogen atom, which is in good agreement 
with values obtained for diffuse ISM in other studies. The 
mass obtained for L1780 is 18M0. 

- The relation between Ay and T(200/im) shows nonlinearity 
at high Ay values, indicating that in the cloud core the far- 
infrared dust emissivity has increased by a factor of ~1.5. 

- The comparison of the 80/im and the lOOyum emission with 
the emission at longer wavelengths indicate that not only 
the 80;um, but also the lOOyum and even the 120jum emission 
may contain some small-grain contribution. 
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